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Abstract 

Temperature dependent Ferromagnetic Resonance measurements performed as a 
function of diameter on Nickel nanowire arrays reveal several interesting features in 
these systems. With diameter decrease from 100 nm to 15 nm, a transition induced 
by surface anisotropy increase is observed at 50 nm in easy axis orientation from 
parallel to perpendicular with respect to individual nanowire geometric axis. Anal- 
ysis of resonance field H res temperature variation (between liquid Helium and room 
temperature) reveals underlying strong magneto-elastic effects in small and large 
diameter nanowire arrays with potential applications in recording and spintronics. 
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1 Introduction 



Ferromagnetic nanowires possess interesting properties that might be exploited 
in current or future spintronic devices such as race-track type magnetic non- 
volatile memory called MRAM (based on transverse domain-wall dynam- 
ics [1,2]) and magnetic logic devices [3,4,5]. They might also be used in partic- 
ular (dissipation-less) magnonic devices [6] extending microwave components 
to the nanoscale regime [7]. 

Mermin- Wagner [8] theorem forbids (Heisenberg-type) magnetism in low- 
dimensional systems (for a dimension < 2) with short-range interactions. 
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Thus, ferromagnetic nanowires being a quasi-one dimensional system display- 
ing magnetic properties represent an interesting system from the fundamental 
point of view. 

In addition to their interest in fundamental magnetism, ferromagnetic 
nanowires have many applications in microwave components such as circula- 
tors [7], superconducting single-photon GHz detectors and counters [9], infor- 
mation storage (as recording media and read- write recording devices), Quan- 
tum transport (Giant magnetoresistance [10]) as well as in Quantum comput- 
ing and Telecommunication devices. 

They are simpler than nanotubes since their physical properties do not de- 
pend on chirality and they can be grown with a variety of methods [11]: 
Molecular Beam Epitaxy, Electrochemical methods (Template synthesis, An- 
odic Alumina filters), Chemical solution techniques (Self-assembly, Sol-Gel, 
emulsions...) and can be grown with a tunable number of monolayers and 
length [12]. 

Ordered arrays of nanowires may be of paramount importance in areas such 
as high-density patterned media information recording an example of which 
is the Quantum Magnetic Disk [13]. They might be also of interest in novel 
high-frequency communication or signal-processing devices based on the ex- 
ploitation of spin-waves (in magnonic crystals made of magnetic superlattices 
or multilayers) [6] to transfer and process information or spin-currents with 
no dissipative Joule effect. 

In this work, Nickel ferromagnetic nanowire arrays (FNA) are fabricated with 
an electrochemical deposition method [14] similar to the one used by Kartopu 
et al. [15]. During deposition, the time dependence of electrical current is 
carefully monitored and recorded since the current intensity profile versus 
time [11] reveals underlying growth mechanisms (coverage of pore walls, filling 
of pore interior by growth, growth at the outbound pore end, hemispherical cap 
growth over nanowire tips, percolative growth outside pores and formation of 
a 2D film) allowing us to finely tune the nanowire growth process and control 
it. 

After structural and magnetic characterization [16], angle dependent FMR 
measurements in the X-band (9.4 GHz) are performed to extract the effective 
anisotropy field H e g versus angle while varying temperature from liquid He- 
lium (4.2 K) to room temperature and changing nanowire diameter from 15 
nm to 100 nm. The length of the Nickel FNA is 6 fim for all diameters and 
the average interwire distance is 350 nm. 

We find that surface anisotropy (for small diameter) and underlying magneto- 
elastic effects play an important role (for all diameter) in those systems that 
might be exploited in novel storage devices, magnonic and spintronic compo- 
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Fig. 1. Magnetization M, applied field H and corresponding angles 0, (j), Oh Ah 
they make with the nanowire axis that can be considered as an ellipsoid-shaped 
single domain with characteristic lengths a = d/2 and c with d the diameter. When 
the aspect ratio c/a is large enough the ellipsoid becomes an infinite cylinder. 



This work is organized as follows: In section 2, FMR measurements are pre- 
sented and later analyzed in section 3. We conclude the work in section 4 and 
Appendix A details the FMR angular fitting procedure, Appendix B covers 
dipolar effects while Appendix C deals with transverse single domain issue. 



2 FMR Resonance field versus field angle and temperature 

Individual wires inside the array are aligned parallel to each other within a 
deviation of a few degrees. They are characterized by a cylindrical shape with 
a typical variation in diameter of less than 5% with a low-surface roughness 
and a typical length of 6 microns. 

FMR experiments are performed with the microwave pumping field h r f oper- 
ating at 9.4 GHz with a DC bias field H making a variable angle 9 H with the 
nanowire axis. 

Previously, several studies have considered reversal modes by domain nucle- 
ation and propagation (see for instance Henry et al. [17] for an extensive 
discussion of the statistical determination of reversal processes and distribu- 
tion functions of domain nucleation and propagation fields). Moreover, Ferre 
et al. [18] and Hertel [19] showed the existence of domains with micromagnetic 
simulations). We do not consider domain nucleation and propagation in this 
work and rather concentrate on transverse single domain case as explained in 



nents. 
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Appendix C. 



Thus, the angular dependence of H res in the uniform mode is obtained by con- 
sidering an ellipsoid with energy E comprised of a small second-order effective 
uniaxial anisotropy [20] contribution K\ and shape demagnetization energy. 
Their sum is the effective anisotropy energy Ea to which we add a Zeeman 
term Ez due to the external field H: 



E = E A + E z = (Ki + ttM|) sin 2 9 

— M S H [sin 9 sin 9 H cos(0 — 4> H ) + cos 9 cos 9 H ] 



(1) 



9 is the angle the magnetization makes with the nanowire axis (see fig. 1). The 
resonance frequency u r is obtained from the Smit-Beljers [21] formula that can 
be derived from the Landau-Lifshitz equation of motion with a damping term 
a. Calculating the angular second derivatives of the total energy: 
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This provides a relationship between the effective anisotropy field -tf e // and 
the external field H at the resonance frequency. 

Theoretically, the effective anisotropy field -tf e // can be obtained from the 
vectorial functional derivative of the energy Ea (eq. 1) with respect to magne- 
tization H e fj = — H|- that becomes in the uniform case the gradient with 
respect to the magnetization components i? e // = — fff-- 

The frequency-field dispersion relation obtained from the Smit-Beljers equa- 
tion is: 
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1 + a 2 ) [H eff cos 29 + H cos(# - 9 H )\ 
x y/[H eff cos 2 9 + H cos(9 - 9 H )} 



(3) 



where H e ff = ^ ss . At the resonance frequency u> = u> r , 9 = 9 H and the 
applied field H = H res in the saturated case. In the unsaturated case the 
magnetization angle 9 ^ 9h and one determines it directly from energy mini- 
mization. 



Calling at equilibrium, magnetization orientation 9q (taking 



bu) we de- 
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termine it by evaluating the derivative (^f ) = 0. We get: 

K ef f sin 20 = M S H sin (0 H - 9 ) (4) 

where K eff = K 1 + irM^. 

Equations 4 and 3 are used simultaneously to determine the resonance field 
H res versus angle H at any temperature (see Appendix A). 

2. 1 General analysis of experimental results 

We add to the previously defined anisotropy energy Ea, the dipolar interaction 
among nanowires yielding i? e // as the sum of the demagnetization field H dem = 
2n Ms, the dipolar interaction field H d i P and the magnetocrystalline anisotropy 
field H K : 

IK i 

Heff = H dem + H dip + H K = 2 7r M s + H dip + — — (5) 

M S 

The dipolar field depends on porosity P (filling factor) in a way such that the 
demagnetization and dipolar fields (see Appendix B) are written as a single 
term 2ttM s (1 -3P). 

Experimentally, the resonance field H res peaks [22] at oo r /l, hence it is possible 
to extract the effective anisotropy field H e ff through the use of eq. 3. 

From the measured resonance field H res versus field angle 6h the (/-factor, sat- 
uration magnetization Ms and cubic anisotropy constant K\ are determined 
with a least-squares fitting method (see Appendix A). 

This yields the following table 1 containing fitting parameters K 1 and M s 
(Anisotropy and saturation magnetization) versus diameter. 

From table 1, one infers that as the diameter increases the Ni bulk values are 
steadily approached which is a good test of the FMR fit. 

2.2 Analysis of FMR results at room temperature 

Measured FMR absorption derivative spectra [23] are similar for (50, 80 and 
100 nm) diameters but differ from the 15 nm case. Angular H res curves versus 
Oh show a minimum at Oh = 0° for the large diameters (50, 80 and 100 nm) 
and a minimum at H = 90° for the 15 nm case (see Fig. 2). 
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Table 1 

Room temperature fitting parameters K\ and Mg with corresponding Nickel 
nanowire diameter d and average separation D. Effective H e ff and anisotropy Hk 
fields are determined with Smit-Beljers. Comparing with bulk Nickel parameters at 
room temperature: K\ = —4.5 x 10 4 erg/cm 3 , K 2 = 2.3 x 10 4 erg/cm 3 , Ms=485 
emu/cm 3 we infer that as the diameter increases we get closer to the bulk values as 
expected with K\ changing by about two orders of magnitude (see Appendix A). 

From this angular variation, we infer that 15 nm samples behave differently 
from larger diameter samples with a transition observed about 50 nm in agree- 
ment with the hysteresis loop VSM measurements (as displayed in fig. 3). 

Previously Nielsch et al. [24] mentioned a change of behavior in nanowire 
arrays because of the existence of a coherence diameter in Ni to be ~ 40 nm. 

The coherent diameter d c separates coherent (Stoner-Wohlfarth style or ho- 
mogeneous) from inhomogeneous reversal (reversal by curling) obtained by 
equating nucleation fields in both cases [25]: d c = 2^~m^ where A is the 
exchange stiffness constant (for Ni, it is about [26] 1.5xl0~ 6 erg/cm). 

Approximating the nanowire by an infinitely long cylinder, the demagnetizing 
factor along the minor axis N a — |, thus: d c = 2y^-^ with q = 1.84 the 

first positive zero of the first kind Bessel function J\{x) derivative dJ ^ ■ Con- 
sequently Nielsch et al. [24] estimate is recovered. Nevertheless in our case, 
the change in behavior is probably due to change in the values of anisotropy 
constant K\ and saturation magnetization Ms with the diameter as seen in 
table 1. 

Surface anisotropy might be responsible for this transition since the FMR fit 
displayed in fig. 1 enables us to define a diameter dependent anisotropy of the 
form K\ = Ky + ^ where the volume anisotropy Ky = —2.89 x 10 5 erg/cm 3 
and the surface anisotropy coefficient Ks = 3.25 x 10 6 erg.nm/cm 3 with the 
diameter d expressed in nm. The 1/d dependence is akin to thin films where 
d is the film thickness [27] . 

The transition has also been observed by Kartopu et al. [15] but left with no 
explanation. Generally, EA orientation of Nickel FNA is believed to be de- 
termined by the interplay of volume and shape anisotropies. This implies two 
regimes exist: for small diameter d (< 50 nm) the EA is along the nanowire 
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Fig. 2. Measured resonance field H res (Ted error bars) versus angle Oh at a frequency 
of 9.4 GHz and at room temperature. The fit with equations 3 and 4 is the continuous 
curve (black line) as shown for all nanowire diameters 15, 50, 80 and 100 nm. The 
angle 9h for which H res is minimal gives the easy axis orientation: it is 90° for 15 
nm and 0° for all larger diameters. 

axis (Isolated nanowire case) whereas for large diameter d (> 100 nm) it is 
perpendicular to it (similarly to the thin film case). 

Nonetheless, a third regime occurs when d is very small (< 35 nm) and sur- 
face anisotropy intervenes with the EA orientation perpendicular anew to the 
nanowire axis. This explains our results for d < 15 nm as well as Kartopu et 
al.'s [15] for <i < 35 nm. 

Considering the fit parameters displayed in table 1 one may understand the 
easy axis direction change as follows. In the 15 nm diameter case, the values 
of K\ and Ms are respectively -1.909 x 10 6 erg/cm 3 and 988.22 emu/cm 3 
yielding the ratio as 1.94 meaning that magnetocrystalline anisotropy is 
more important than shape anisotropy that tends to align the magnetization 
along the nanowire axis. Moreover, since K\ < 0, the anisotropy is planar and 
consequently the easy axis is in the plane perpendicular to the nanowire axis. 
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Fig. 3. (Color on-line) Room temperature VSM measured hysteresis loops M/M$ 
versus H (in Oersteds). The field H is perpendicular (blue for _L) or parallel (red 
for / /) to the wire axis in the 15, 50, 80 and 100 nm diameter cases. The transition 
of easy axis orientation from perpendicular (_L) in the 15 nm case to parallel (/ /) 
for diameter > 50 nm is clearly visible in the loops (loop merging at 50 nm and 
interchange for larger diameters) and confirms what is observed with FMR in fig. 2. 



In the 50 nm case, the corresponding values are respectively -1.621 x 10 5 
erg/cm 3 and 451.95 emu/cm 3 yielding a ratio of 0.79. This means mag- 
netocrystalline anisotropy is less important than shape anisotropy that tends 
to align the magnetization along the nanowire axis. 

In fig. 4 the theoretical behavior at room temperature of the resonance fre- 
quency tc> r /7 versus field is shown. The intersection of the FMR measurement 
frequency line at 9.4 GHz we infer that in both diameter cases, we have a low 
field mode for 9h = 0° and a high field mode at 9h = 90° as observed in fig. 2. 
In the latter the low-field mode is around 1000 Oe for the 15 nm diameter and 
1500 Oe for the 100 nm case. In contrast, fig. 4 indicates smaller values (250 
Oe for the 15 nm diameter and 500 Oe for the 100 nm case) stemming from 
averaging effects arising from a dispersion of array anisotropy or geometry. 
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Fig. 4. (Color on-line) Room temperature theoretical resonance frequency ^ as 
a function of field for the 15 nm (left) and 100 nm (right) diameter with angles: 
6h = 0° and Oh = 90° . The horizontal line is the FMR measurement frequency 
of 9.4 GHz. The low-field quarter of a circle shaped curve is the unsaturated case 
where the equilibrium magnetization angle is different from Oh- It can be expressed 
as ^ = l H^ff ~ H 2 . In the 15 nm, we have a high field mode for 9jj = 0° and a 
low field mode at Oh = 90° as shown in fig. 2. It is exactly the opposite for diameter 
d >15 nm as observed in fig. 2. 

2.3 Analysis of FMR results versus temperature 



Extending the above approach to the temperature dependent case by explicitly 
expressing Ms and K\ as M$(T) and K\(T) (taken as the first cubic anisotropy 
constant [28]), such that the temperature dependent effective field writes: 



H eff (T) = 2 7rM s (T)(l-3P) 



2Ki(T) 
M S {T) 



(6) 



The temperature variation of the resonance field H res (T) is obtained from the 
Smit-Beijers equation 3 (neglecting a) after relating it to the effective field 



H. 



eff- 



When H is applied along or perpendicularly to the wire axis, we get: 
H res (T) 



J + H eff (T) 



6 H = 0° 



H res (T) = - 









H eff( T ) ~ \ 




. 7 . 



9 H = 90° 



(7) 



Performing the analysis of measured FMR lineshapes versus field, tempera- 
ture and angle (fig. 5 displays the FMR lineshapes as a function of field for 
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Fig. 5. (Color on-line) Measured FMR absorption lineshapes as a function of field 
for temperature varying from 4.2 K to room temperature for the 100 nm nanowire 
array. The field angle is 6 H = 0° (top) and H = 90° (bottom). 

temperature varying from 4.2 K up to 300 K for two field angles Oh = 0° and 
9h = 90° in the 100 nm case) entails extracting H res (T) (as displayed in fig. 6) 
and analyze its behavior with theoretical models. 



In Fig. 6 theory and experiment of the temperature dependence of H res (T) are 
displayed. The disagreement observed indicates that additional temperature 
dependent anisotropics might intervene in the behavior of H res (T). 

Noticing that in all cases H res (T) decreases with temperature a possible source 
of discrepancy might originate from the difference between thermal expansion 
coefficients of the metallic nanowire am with respect to the supporting di- 
electric material (DM) a DM . Thus, magneto-elastic effects might play an im- 
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Fig. 6. (Color on-line) Experimental (crosses) and theoretical (continuous lines) 
resonance field H res as a function of T £ [<i.2K, 300K] for the 15 nm (left) and 
100 nm (right) diameter with both angles: Oh = 0° (red crosses for experiment and 
green lines for theory) and Oh = 90° (blue crosses for experiment and magenta lines 
for theory). Magneto-elastic terms are not considered in anisotropy. Larger values 
of H res occur for Oh = 0° in the 15 nm case and for Oh = 90° in the 100 nm case. 



portant role in these systems. Moreover, the difference between am and a^M 
might be strengthened further with chemical bonding or interaction effects 
between nanowires and the supporting material. 

In order to account for the effect of temperature effect, let us start from 
the definition of effective anisotropy constant comprising magnetocrystalline, 
shape and dipolar contributions K e ff — K\ + 7rM|(l — 3P) and extend it in 
a way such that it becomes function of temperature. 

Taking account of surface and magneto-elastic contributions, yields: 

K eff (T) = K X (T) + ttM|(T)(1 - 3P) + K me (T) + (8) 



K\(T) and K me (T) are respectively the magnetocrystalline [29] and magne- 
toelastic [30]. In the following, temperature dependence for all anisotropies is 
determined except K$- 

Regarding saturation magnetization M$(T) we use Kuzmin [31] parameteri- 
zation, i.e.: 

M S (T) = M 5 (0)[1 - sx 3/2 - (1 - s)x p } 1/3 (9) 



with x = T/T c , Ms(0) = 57.6emu/g, p = 5/2, s = 0.15 and the Curie temper- 
ature [31] T c = 628K. 
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The temperature dependence of the magnetocrystalline coefficient K\(T) fol- 
lows the Carr model [28], i.e.: 

10 

(10) 

with 5 = 1.74 and K = —4.5 x 10 5 erg/cm 3 Ni first anisotropy constant at 
T = OK which is increased by one order of magnitude with respect to its room 
temperature [28] value. 

Magnetoelastic anisotropy K mc (T) is estimated along the following: K mc (T) = 
|As<r with a = E Y e(T) where a is stress, E Y is Young modulus and e(T) the 
temperature dependent deformation. 

Deformation e(T) is estimated from the difference of dilation coefficients be- 
tween Nickel and the DM: e{T) = (a m - a DM )AT with AT the difference 
between fabrication (300 K) and measurement temperatures, am, «dm are 
respectively the thermal expansion coefficients of Nickel and of the DM. Their 
actual values are: a Ni = 13. x 10~ 6 /K and in the case Alumina (A1 2 3 ) is the 
DM, a DM = 8.1 x 10~ 6 /K. 

Xs is Nickel saturation isotropic magnetostriction constant [32] equal to —34.10 
Young modulus of bulk Nickel is about 2.05 x 10 12 dyn.cm~ 2 yielding a tem- 
perature dependent magnetoelastic anisotropy varying as: K mc (T) ^512 AT. 
This result being approximate, we use an expression of the form K mc (T) « 
a AT and searched for the best value of a that might explain all results for 
all diameters and field angles. We found the best value as 1500 erg/K.cm 3 as 
displayed in fig. 7. 

Theoretical and experimental values of the resonance field H res versus temper- 
ature are compared in fig. 7 after accounting for the magneto-elastic anisotropy 
K me (T) contribution in the form K me (T) 1500 AT. The agreement as a func- 
tion of temperature for both nanowire diameters and both angles {6 H = 0° and 
9jj = 90° ) of the field is a strong indication of the presence of magneto-elastic 
effects. 



3 Discussion and Conclusion 

We have performed angle and temperature dependent FMR on Ni nanowire 
arrays with variable diameter and shown with FMR that the easy axis orienta- 
tion for the 15 nm diameter sample is perpendicular to the wire axis in sharp 
contrast with the 50 nm, 80 nm and 100 nm samples. Note that we expect 



K\(T) — Kq(1 — 8x) 
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Fig. 7. (Color on-line) Experimental and theoretical resonance field H res accounting 
for magneto-elastic contribution as a function of T £ [4.2K, 300K] for the 15 nm 
(left) and 100 nm (right) diameter with both angles: 9h = 0° (red crosses for ex- 
periment and green lines for theory) and 6h = 90° (blue crosses for experiment and 
magenta lines for theory). The magneto-elastic contribution improves substantially 
agreement between theory and experiment. 



(from bulk Ni) the easy axis along the wire axis by comparing the value of 
shape energy with respect to anisotropy energy. 

Results obtained from the angular behavior of H res versus 6h show that H res 
is minimum at 90° for the 15 nm sample whereas it is minimum at 0° for 
the larger diameter samples agree with hysteresis loops obtained from VSM 
measurements and confirm presence of the transition of easy axis direction 
from perpendicular at 15 nm to parallel to nanowire axis at 50 nm diameter. 

The transition observed at 50 nm is interesting because of several potential 
applications in race-track MRAM devices. Yan et al. [4] predicted that in 
Permalloy nanowires of 50 nm and less, moving zero-mass domain walls may 
attain a velocity of several 100 m/s beating Walker limit obeyed in Permalloy 
strips with same lateral size. Hence, nanowire cylindrical geometry in contrast 
to prismatic geometry of stripes bears important consequences on current 
injection in nanowires that applies Slonczewski type torques [33] on magneti- 
zation affecting domain wall motion with reduced Ohmic losses [34] . 

At low temperature, we find that H e ff increases when temperature is de- 
creased. This may be attributed to the increase of Ms(T), as temperature is 
decreased, affecting all anisotropy fields (H K , H dem , and H dip ) that depend on 
M S (T). 

Separating the various contributions to anisotropy by using thermal, fre- 
quency, diameter and angular variations we have been able to pinpoint the 
main contribution to anisotropy at low temperature as stemming from magneto- 
elastic effects between the nanowire array and the DM supporting it. 
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Magneto-elastic effects affect directly the thermal variation of H e ff that in- 
crease at low temperature due to the increasing difference between the thermal 
expansion of the metallic nanowire and the DM. Thermal effects must be an- 
alyzed properly for building extremely high density storage devices. 

Ordered arrays of nanowires are good candidates for patterned media and 
may also be used in plasmonic applications such as nano-antenna arrays or 
nanophotonic waveguides in integrated optics [35]. Recently [36], heat assisted 
magnetic perpendicular recording using plasmonic aperture nano-antenna has 
been tested on patterned media in order to process large storage densities 
starting at 1 Tbits/in 2 and scalable up to 100 Tbits/in 2 . 

In thermally-assisted perpendicular magnetic recording, a waveguide delivers 
light to a plasmonic nano-antenna placed just above the disk platter surface 
creating an intense optical pattern in the near-field region. However, this heats 
the disk on the nanometer scale (around 25 nm) and consequently tempera- 
ture dependent magneto-elastic effects analyzed in this work might affect the 
physics of high-density information writing in these systems. 
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4 Appendix A: Angular FMR fitting procedure 

We have developed a procedure based on a least squares minimization proce- 
dure of the curve H res versus 8h to the set of n experimental measurements 
[xi, Ui\ i=ln where Xi = 9 Hi and = H res ((3; xt). (3 represents the set of 
parameters to fit M s , K\. 

Hence the set of minima equations for the data points are: 

Hi] 2 minimum, 

0°)| minimum, (11) 
90°) | minimum 

where the values y Q (resp. yg ) are the experimental values of H res correspond- 
ing to the angle 6 H = 0° (resp. 9 H = 90°). 



-Y}H res ((3;Xi) - 
n i=i 

1 2/o — H res {f3; Xi = 
2/90 - H res ((3;xi = 
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The fitting method is based on the Broyden algorithm, a generalization to 
higher dimension of the one-dimensional secant method [37] that allows us to 
determine in a least-squares fashion, the set of unknowns g, M s ,Ki. Broyden 
method is selected because it can handle over or under-determined numerical 
problems and that it works from a singular value decomposition point of view 
[37]. This means it is able to circumvent singularities and deliver a practical 
solution to the problem at hand as an optimal set [37] within a minimal distance 
from the real one. A practical solution might artificially increase some of the 
fitting values as observed in table I. 



5 Appendix B: Dipolar interaction field 

Demagnetization energy is given [38] by E dem = 2irN i jM i Mj where the de- 
magnetization coefficients are such that < < 1 and i,j = 1,2,3. Since 
the demagnetization field is given by H dem = — 9 Q^p , component % of H dem 
is —AirNijMj. We have two limits: 

• Single isolated nanowire (for which P = and N xx = 1/2, N yy = 1/2, N zz = 
0), E dem = ir{M* + Ml) = 7r(M| - M 2 Z ) and H dem = 2nM s cos 9z with z 
the unit vector along z. 

• Thin film limit (for which P = 1 and N xx = 0, N yy = 0, N zz = 1), E dem = 
2nMl and H dem = -4nM s cos 6z. 

When the sample is saturated along z (nanowire axis) M x = 0, M y = 0, M z = 
Ms, the demagnetization field in both cases has a single z component: H dem = 
2-jiMs in the single wire case and H dem = —47rMs in the thin film case. Making 
a linear interpolation between these two limits, we get: H dem = 2nM s (l — 3P). 

If we rather consider a 2D square lattice of nanowires with parameter D the 
average nanowire separation, the porosity is given by: P = ^ and we can 
calculate directly the dipolar energy as follows. 

Starting from a single dipole po surrounded by an array of dipoles on the 
lattice pi, the interaction energy is: 

p _^ PoPi 3(p -r)(pi-r) 

&int - 2^ 3 5 
i 

This is equivalent to a field H dem acting on the dipole such that: E int = 
— po ■ H dem . 

When po is directed along the z axis, p = (0, 0, 1) and p • H dem = H dem z 
with H dem , tZ the z component of H dem . 
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(12) 



At any node of the lattice Pi = (p x ,Py,Pz), and r = (iD,jD, 0), hence 
Po-r = 0,pi-p = p z . 

The field H dem ,z is then: 

Hdem ' z = ~ . ? ^ = ~ . ? n3f ? -2 I ^1 ( 13 ^ 

i,j=— oo !j=-oo V 6 'J J 



When all surrounding dipoles are saturated along z, p z = MsV with V = na 2 l 
the nanowire volume : 

M S V +2? 1 A.2M S V . I . 

= JL(?T?jI K = - 4 ' 2MsP( 5' (14 > 



This should be compared to the interpolated dipolar term — QttM s P. 

The evaluation of Hd em , x the x component of Hd em proceeds along the same 

lines. 

Taking p along the x axis, p = (1,0, 0),p = (p x ,p y ,p z ) and r = (iD,jD,0), 
the dipolar energy writes: 

F n H H V ^p x -3D^p x +yp y ) 

J^int — —PO ■ -n-dem — —J^dem,x — ~ 2^ 5 v 10 / 

ij=-oo r 



Therefore: 



rr Px(f±f) -3(Pp x +ijp y ) 

ndem, X — 2^ m /.n , - 2 \5 \ ' 

iJ=-00 D*(l 2 +f)2 



If all the surrounding dipoles are saturated along x,p x = M s V,p y = p z = 0, 
we get: 

Hlmx =^L i ^.jjjw: (17) 



Thus we obtain: 

if^s = 2.\M S P{^) = - l -Ede m ^ (18) 
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6 APPENDIX C: Single domain radius 



A simple argument given in Landau-Lifshitz Electrodynamics of Continuous 
Media [38] gives the size below which a single domain (or Stoner-Wohlfarth) 
behavior is observed in the plane perpendicular to nanowire axis. It is based 
on the following: When the demagnetization energy 2irN i jM i Mj ~ 2nN c Mg 
(where N c is the demagnetization coefficient along some preferred axis, usually 
the long one in an ellipsoid- approximated nanowire) is equal to the exchange 



R. 



(A^, k — 1, 2, 3 is the exchange stiffness constant 



- A a typical exchange stiffness constant (regardless of 
■ Exchange energy is the largest contribution to 



along i,j directions). 
Considering that Aij 
results in R S( i ~ 

non-uniformity energy due to spatial variation of the magnetization M. 
When the change in the direction of M occurs over distances that are large 
compared to interatomic distances, non-uniformity energy can be expressed 
through derivatives of M with respect to spatial coordinates (see Landau- 
Lifshitz [38] and Brown [39]). Exchange stiffness constant A^ is on the order 
of Heisenberg exchange energy per unit length J / ao (ao is the average nearest 
neighbour distance in the nanowire material). 

10000 



woo 



100 




w 



Single Domain 
Coherent 
Landau 



0.05 0.1 



0.15 0.2 



0.25 0.3 0.35 



Fig. 8. Single domain radius (nm) for a Nickel prolate (elongated) ellipsoid as a 
function of the demagnetization coefficient along its axis N c . The Landau approxi- 
mation R S( i ~ \J 2-kN nT 1 ls sca l e d i n a wa y such that it agrees with the exact result 
when m — > 0. The coherent radius, in contrast, increases with iV c and is around 
10 nm. The room-temperature Nickel data taken from Kittel [40]) are M s = 485 
Gauss; ao = 2.49 A. The exchange stiffness constant A is taken as 10~ 6 erg/cm. 

Typically J ~ 10 meV and ao ~ 1A, hence we get A^ ~ 10~ 6 erg/cm (see 
fig. 8). The exchange length £ ex is defined as £ ex = \J ' A/K with A, the exchange 
stiffness constant (A ~ 10~ 6 erg/cm). It is obtained from the comparison 
between the exchange energy j^^r^r ~ and the anisotropy energy 
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^MiMj ~ K with K representing the anisotropy constant of the material. 
In soft materials whose anisotropy constant K ~ 0, one uses the magnetostatic 



exchange length defined as (i ex = ^A/M^ since K ~ in soft materials. It 
is obtained from the comparison between the exchange energy and the 
demagnetization energy 2tiN c M 2 . 

The length i ex = \J~A/~K is in fact on the order of the domain wall thickness, 
therefore we follow Frei et al. [41] to estimate rigorously the single-domain 
radius R s d from a minimization of the energy using Euler variational equations. 

The energy density of an infinite cylinder (see fig.l) made of a ferromagnetic 
material accounting for only exchange terms and Zeeman energy is given in 
cylindrical coordinates (r, if, z) by: 

,JM, 2 f dM 2 dM 2 ,^ „ 

^~dr^ + Vdip + ~dz~ ~ M.H (19) 



The total energy is found from an average over the volume: 

1 J r . r ,d6 2 sin 2 ^ 
' J 

«o/2 



i ,. rlf) 2 sin 2 f) 

E = — 2 {A[(-) + ]-M s Hcos(6- ( f ) )}2Trrdr (20) 
7Tcr J dr r z 



Minimization of the energy obtained from nulling the variational derivative 
with respect to the angle 9 leads to: 

d 2 6 1 d6 1 

- 2— - 2-— + 2- sin(0) cos(^) + M S H sin(0 - 0) = (21) 
dr l r dr r 



When we consider small values of 9, the above differential equation becomes 
a Bessel equation in terms of 9(r). 

In order to evaluate the average of the different energy terms, we replace 6{r) 
by 7r/2 in the evaluation of the exchange energy term that becomes Ex = ^r- 
This can be simply averaged over the ellipsoid volume (with short axis 2a and 
long axis 2c see fig. 1) representing the nanowire: 



r 2 

S =i^/ Wr / d 4 (22 > 

3 a /2 



a is a cutoff length given by the nearest-neighbour distance or the lattice 
parameter of the corresponding crystal (see fig 8). The average magnetostatic 
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(demagnetization) energy can be calculated following ref [42] accounting for 
the curling of the magnetization from the initial state with M : 



with N c (the long axis demagnetization coefficient). 

Comparing the average exchange and demagnetization energies gives: 



Solving eq. 24 for R sd , in the case of Ni, fig. 8 gives the variation of R s d versus 
N c . From the figure, we infer that for Nickel nanowires R sc i is within a few 100 
nm range. 
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